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422. 


ON THE GEODESIC LINES ON AN OBLATE SPHEROID. 


[From the Philosophical Magazine, vol. xL. (1870), pp. 329—340.] 


THE theory of the geodesic lines on an oblate. spheroid of any excentricity what- 
ever was investigated by Legendre(*?); and the general course of them is well known, 
viz. each geodesic line undulates between two parallels equidistant from the equator 
(being thus either a closed curve, or a curve of indefinite length, according to the 
distance between the two parallels): at a point of contact with the paraliel the curve 
is, of course, at right angles to the meridian; say this is V, a vertex of the geodesic 


Fic. 1, 
C 


line, and let the meridian through V meet the equator in A; the geodesic line proceeds 
from V to meet the equator in a point N, the node, where AN is at most = 90°; 
and the undulations are obtained by the repetition of this portion VN of the geodesic 
line alternately on each side of the equator and of the meridian. 


1 Mém. de UInst. 1806; see also the Ewer, de Calcul Intégral, t. 1. (1811), p. 178, and the Traité des 
Fonctions Elliptiques, t. 1. (1825), p. 360, 
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I consider in the present paper the series of geodesic lines which cut at right 
angles a given meridian AC, or, say, a series of geodesic normals. It may be remarked 
that as V passes from the position A on the equator to the pole C, the angular 
distance AN increases from a certain determinate value (equal, as will appear, to 
© 90°, if C, A are the polar and equatorial axes respectively) up to the value 90°; 
and it thus appears that, attending only to their course after they first meet the 
equator, the geodesic normals have an envelope resembling in its general appearance 
the evolute of an ellipse (see fig. 1 and also fig. 2), the centre hereof being the point 


B at the distance BA =90°, and the axes coinciding in direction with the equator 
BA and meridian BC: this is in fact a real geodesic evolute of the meridian CA. 
The point a is, it is clear, the intersection of the equator by the geodesic line for 
which V is consecutive to the point A (so that 2 BOA = (1 — A 90°); and the point 
y is the intersection of the meridian CB by the geodesic line for which V is con- 
secutive to the point C; and its position will be in this way presently determined. 
I was anxious, with a view to the construction of a drawing and a model, to obtain 
some numerical results in relation to a spheroid of considerable excentricity, and I 


selected that for which £z} (polar axis =4 equatorial). 


Before proceeding further, I remark that Legendre’s expression “reduced latitude ” 
is used in what is not, I think, the ordinary sense; and I propose to substitute the 


Fie. 3. 
0’ 


0 AE ae 


term “parametric latitude”: viz., in fig. 3, referring the point P on the ellipse by means 
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of the ordinate MPQ to a point Q on the circle, radius OK (= OA, fig. 1), and drawing 
the normal PT, then we have for the point P the three latitudes, 


A. =ZPTK, normal latitude, 
A” =Z POK, central latitude, 
A =ZQOK, parametric latitude; 


viz. A” is the parameter most convenient for. the expression of the values of the 
coordinates s, y (e=Acosr’, y=Csin2’) of a point P on the ellipse. The relations 
between the three latitudes are 


Ld C / C? 
tan X ar | tan À = 4; tan N, 


so that A”, ~, à are in the order of. increasing magnitude. I use in like manner 
l, V, Ù’ in regard to the vertex V. The course of a geodesic line is determined by the 
equation 

cos X sin a = const., 


where 2’ is the reduced latitude of any point P on the geodesic line, and a is at 
this point the azimuth of the geodesic line, or its inclination to the meridian. Hence, 
if V be the parametric latitude of the vertex V, the equation is 


cos X’ sin a = cos l’ 


(whence also, when \’=0, a =90°—l'; that is, the geodesic line cuts the equator at 
an angle =U’, the parametric latitude of the vertex). The equation in question, 
cos X sin a=cos U’, leads at once to Legendre’s other equations: viz. taking, as above, A, C for 


C2 
and to determine the position of P on the meridian, using (instead of the parametric 
latitude X’) the angle ¢ determined by the equation 


the equatorial and polar semiaxes respectively, and 6 for the excentricity, = J, 


sin ’ 
sspe TH 
and writing, moreover, s to denote the geodesic distance VP, and A to denote the 


longitude of P measured from the meridian CA which passes through the vertex V, 
these are 


ds = do VC? + A%& sin? l cos? ¢, 
cos l’ dé VO: + AS sin? l cos? d 


sci A 1 —sin*/’ cos? Š 


which differential expressions are to be integrated from ¢=0; and the equations then 
determine 2’, s and A, all in terms of the angle ¢,—that is, virtually s and A, the 
length and longitude, in terms of the parametric latitude X. 


Cc. VIL. 3 
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Writing, with Legendre, 
A? sin? Y 


C$ @y De aml’ ramet, 
C? . 
k=l- e, = 024 AF sinl’ =1-—&sin*l; 
also 
n=tan?l’, M C q4 


~ A cost’ A cosl’ 
then the formule become 


d= Ç dg v1 ean 6, 


dp V1- esing 


e oian nr TT 


Hence integrating from ¢=0, and using the notations F, E, II of elliptic functions, 
we have 


s= FEC $), 
M 
A =— {(n+e) (n, c, $)—eF(, $); 


viz. these belong to any point P whatever on the geodesic line, parametric latitude 
of vertex =l; and if we write herein ¢=90°, then they will refer to the node N, 
or point of intersection with the equator. 


The position of the point a is at once obtained by writing /'=0: viz. this gives 


c=0, b=1, M= n n=0: the differential expressions are ds=Cd¢, dA = ~ dg. Or 


integrating from ¢=0 to ¢=47, we have s=A. Cir, A=. 4m, agreeing with each 


other, and giving longitude of a = 4 47; or, what is the same thing, 2a0OB=47 (1 - 5). 


Writing in the formule l= 90°, we have c=6, b= A a Zo; whence dA =0, or 
A =const., =47, since the geodesic line here coincides with the meridian CB; and 
moreover s=AH(6, p); viz. this is merely the expression of the distance from © of 


a point P on the meridian CB. But we do not thus obtain the position of the point y. 
To find it we must consider a position of V consecutive to O, say, l’=42—-—¢, where 


e is indefinitely small; n is thus indefinitely large, and the integral II (n, c, $) is not 
conveniently dealt with. But it may be replaced by an expression depending on 
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2 
I (E, g, $), where ~ is indefinitely small; viz. (Legendre, Fonct. Ellipt. vol. 1 p. 69) 


we have 


a 1, Vatan e 
II (n, c, ġ)=F (c, $t g erry le C, $), 


VI — @sin? 
where 


a=(1+n) (1 + 2) 
We thus have 


M are ee è 
= T fare, AAi tan aeg CHAG C, #), 


where, Z being small, 


C = ERIO Pe... Se 
II =. C, $) Naan oE 


-£ —-— " gin? $) do 


m cra =(1- 1) FG, $) + EC, $). 


And expanding also the tan term, we thus have 


Vavetn[ V1—c*sin® $ Vn 
“RENO O+ Tee ET T 


-(e+n)|(1- a +580 A|}, 


2 24, 
-Hoire = (142) pr ENEE ae tatovi aana], 
which, in the term in { } neglecting negative powers of n, becomes 


M ( ,- pa a siate 
Aw on fva $r+bF(c, ġ)— E (c, ġ)— cot o VI -— e sin? $} . 
7 M4 . C o , 1 M 
We may moreover write c=6, b= A? ¢=90°-2X, n= 3 M=e, and therefore = = & 80 
that the formula is 


Aze f . $m +F (c, 90°- N) — E (c, 90° — X) — tan XN V1 -= ¢ cosa} à 


= $m —e {tand’ V1 =e cos N + E (e, 90° — N) — bF (c, 90° —X}, 


i ; KaU : 
where I retain c, b as standing for 4/1 — w1 respectively. 
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Writing herein \’=0, we have 
A =Ẹġr —e(E, c- bF, c), 
where the coefficient E, c— bF, c is 
in th ra is be 2 
0 V1—csin? 6 o V1 — esin? 6 


consequently positive; that is, A, the longitude of the node, is iess than 90°, as it 
should be. Hence in order that A may be =90°, we must have 2’ negative, say, 
A =— yp’, where yp’ is positive; and, observing that we may under the signs Æ, F 
write 90° — x’ instead of 90°+ w’, we thus have 


ho =hor+e{Vl1—Ccos*y tan p’ —E(c, 90° — p) +F (c, 90°-p’)}; 
that is, we must have 
tan p’ V1 — ecos? p = E (c, 90°—p’) — F (c, 90° — u’); 


viz. w is here the parametric latitude (south) of the intersection of the meridian CB 
with the consecutive geodesic line—that is, of the point y. As w increases from 
0 to 90°, the left-hand side increases from 0 to œ; and the right-hand side, beginning 
from a positive value and either attaining a maximum or not, ultimately decreases 
to 0; there is consequently a real root, which is easily found by trial. 


Thus z =4, C=473 (angle of modulus = 60°), b=4; or the equation is 
tan w’ V1 — $ cos yw’ = E (90° — p’) — 4F (90° — pv’). 
Using Legendre’s Table IX., we have 


pl 90° — w E E~4F. tan uw’ 1—3 cos? w. 
0° 90° 1-21105 2°15651 6719 ‘0 

10 80 1°12248 1°81252 6693 

20 70 1:02663 1:49441 6530 

30 60 -91839 1:21253 6153 *3819 

40 | 50 ‘79538 ‘96465 "5542 6278 


so that we see the required value is between 30° and 40°; and a rough interpolation 
gives the value w =37°40. But repeating the calculation with the values 37° and 
38°, we have 


F. E-}F. tan p’ V1- ¥ cos? p. 


‘833879 | 1:035870 ‘57419 
821197 | 1:011849 56823 


whence, interpolating, w = 37° 55’. 
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The semiaxes of the geodesic evolute, measured according to their longitude and 
parametric latitude respectively, are thus Ba, long. of a=45°; By, param. lat. = 37° 55’. 
But measuring them according to their geodesic distance, the equatorial radius A being 
taken =1, we have 


Ba=}łr = ‘78540, 


By = (5 i 1) (E, — E (52° 5’)} = 121106 — 82225 = 38881. 


Reverting to the general formule for s, A, but writing therein A = 1, and therefore 


C=V1—8; writing also = 90° (that is, making the formule to refer to the node 
N of the geodesic line), we have 


vVl]-=8 
po e E 
VI — ẹsin? ? 0% 
V1 —® 
= os] (e+e) (n, €)— CF, e}; 


but for the calculation of the second of these formulae by means of Legendre’s Tables 
it is necessary to express II,(n, c) in terms of the functions Æ, F. 
The proper formula is given in Fonct. Ellipt. vol. 1. p. 187; viz. this is 


eo (n o)=}r +728 A, 0) F,c+ FoF (b, 0)-F,c EQ, 0)-E,c F0, 0), 


where A (b, 0)=V1—Db*sin?@. @ is an angle given by the equation cot@=Vn; we 
have n=tan*/’; consequently @=90°—U’. Substituting this value, except that for 
shortness I retain Æ (b; 0), F(b, 0) in place of E (b, 90°- 7), F(b, 90° — V), we have 

A (b, 0) =V1—2 cos? V’, 


= ¥71—(1— &sin®/) cos? l , =sin l; 


and thence 
tan 0A (b, 0) = cot lsin l= 5 . 
whence 
E 0) mrt lr + Fc F= e+ Fb, -EQ o|- E.cF(b, ay} . 
But 
n+c?=tan?l’ + ẹsin? l = sin? / sec? l. 
Hence 


(n + e) II, (n, c)— eF, c = sin? l {sec° TI, (n, ¢)— &F, c} ; 
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and multiplying this by 
V1 — & _N1 — & cos l 
meosl’ — tarl cosl’ 
the exterior factor is 
VI — & cos ltan? 1 cos l 


tan? 7 RS I A 


and we have 


an paan (sec* UII, (n, c) — &F' ch, 


which is the formula I used in the calculations. It would, however, have been better 
to reduce a step further; viz. we have 


al tan Ų 
sect CTL, nah r cos | t 
V1 —& cos 1 
=i |ir+F,c F +P, 0) -EQ @)|- Ec FQ, o), 
eerk (hr +Fc[F(b, 0)-E(b, 0]—-E,cF(b, 0}+ Po, 
and thence 
NEE, 


sec? l'II, (n, c)— èF c= ON. Cs Lie +F c[V1—®&cosl+ F(b, 0)—E(b, 0))— Eck (b, 6)} ; 


or, finally, 
A=ġr + F, c F(b, 0)-—F ck (b, 0)-E,cF(b, 0)+V1—®&coslF,c. 


It is easy with this expression of A to obtain the results already found for the 
extreme values l’= 0°, V = 90°. 


As Legendre’s Tables have for argument, not the modulus c, but the angle of the 
modulus, say x (that is, siny=c=Ssin/), it is convenient to replace V1—& sin? by 
its value cosy; and the formule thus are 


A=4r+F c [V1 — & cos l+ F(b, 0)—E (b, 0)|—E,cF(b, 0), 


where 


C=siny=Ssinl, tanl’/=V1—Stanl, @=90°-/; 


and in the case intended to be numerically discussed, §=443, V1 —&=}4. I take V” 
as the argument, giving it the values 0°, 10°,... 90°, and perform the calculation as 
shown in the Table. 
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F (b, 0), nat. 
E (b, 0), nat. 


47151 | 20548 18686 | 1-6050 


| 
59°37 | -40425 | 22814 | 16532] 1-6910 


2:08962 | 1:03762 
1-48840 | 1-02962 
49°10 | -32737 | 25478 | 14167] 1-7980 1:16024 | -95214 
"92141 | 82827 
*71820 | 67903 


-53083 | -51655 


| 


41°93 | 20989 27626 | 12163] 1-8891 
37:02 19349 29935 | 10645 | 1-9923 
33°68 | 13866 | 31479 | 09557 | 2:0644 


31°57 | -08954 | 32540 | 08850 | 2°1154 35099 | *34716 


17475 | 17431 


_ 

© 

@ 

© 

pa 

So 

tbo 

on 

-_ 

orks) 

-~J 

or 

-~J ! 
wo 

bo 

or 


30°38 | -04387 | 33169 | 08446 | 2-1463 


08316 


* * 

ee t ; q a 

-g S Pa S ‘ . E 

s| 8 | 48 a ejes j 

i a 3 Sf ; Se) & WE A [yy] 2 l 

rah ae: oe e | 3 á |53 $ 

3% a | a Z i 3 

| 5 ~AN 

| o o ve 4 

0 l 45 -7854 
10152308 | -18272 | -38820 2-4446 | -32007| -50693 | 3-2132 45 58 1:6058 | -8029 
20} -86318 | 193610 | -16424 1-4596 | -17272 | 33804 | 2-1779 48 51 1:7012 | -8506 
30] 53547 | 172874 198352 -9627 | -06455 | 20622 | 1-6078 53 2 1-8333 | -9166 
40] -34903 | 1-54286 | 1-81912 -6594 f 1-96445 | -08608 | 1-2192 57 56 1:9801 | -9900 


50} -23266 | 1-36672 | 1-66607 -4635 | 1-85625 | 1-96270 | -9177 
60} 15292 | 1-18446 | 1-49925 | -3157 | 1-72496 | 1-82053 | -6615 
70} 09327 | 2-96974 7-29514 | ‘1973 | 154529 | 1-63379 | +4303 
80] -04431 | 2-64650 | 297819 | -0951 | 1:24242 | 1-32688| +2123 


x, 90°—x in degrees and decimals of a degree, to correspond with Legendre’s Tables, 


Where the columns marked with an * show respectively the longitude of the node, 


and the length (or distance of node from vertex), for the geodesic lines belonging to 
the different values of the argument V’. 
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The remarks which follow have reference to the stereographic projection of the 
figure on the plane of the equator, the centre of projection being the pole (say the 
South Pole) of the spheroid. It is to be remarked that if a point P of the spheroid 
is projected as above, by means of an ordinate into the point Q of the sphere radius 
OK (= OA), then projecting stereographically as to the spheroid and the sphere from 
the south poles thereof respectively, the points P and Q have the same projection. 
And it is hence easy to show that an azimuth a at a point of the meridian 
(parametric latitude M, normal latitude à, and therefore tan N = ps tan r) is projected 
into an angle (a) such that 


sin 0’ 
tan @)n aaa tan a. 

In fact in fig. 3, if we take therein OK, OC for the axes of æ, z respectively, 
and the axis of y at right angles to the plane of the paper, and if we have at P 
on the surface of the spheroid an element of length PR at the inclination a to the 
meridian PK, then if æ, y, z are the coordinates of P, and «+ ôs, y+ òy, 2+ 8z those 
of R, we have 


dv= pcosasin À, 
dz = — p COS 4 COS À, 
y= psina, 
and thence 
tan a= W. Bol 
V 8a? + 82? 


Now, if the meridian and the points P, R are referred by lines parallel to Oz to the 
surface of the sphere radius OA, the only difference is that the ordinates z are 
increased in the ratio C : A; so that if the projected angle be (a), we have 


and then projecting the sphere stereographically from its south pole, the angle in the 
projection is = (a«a). And according to the foregoing remark, the angle (a) thus obtained 
is also the projection of a from the south pole of the spheroid. We have thus 


tan (a) _ 8a? + 82 _ Msin? r+ cos? _. f/l+cot?X —_ sind’ 
tana Gaiti {3 ] 2 f rano V boot aT] A 
Ne +5 bet a/ sint 0+ S costa 


which is the required relation. 


The foregoing equations, 


cos’ sin a = cos l, tanN = $ tan A, 


+ 


tan (a) = = 3 tan a, 


A 
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determine in the stereographic projection the inclination (a) to the radius, or projection 
of the meridian, of the geodesic line (parametric latitude of vertex =/’) at the point 
the parametric latitude of which is =)’; viz. they enable the construction (in the 
projection) of the direction of the successive elements of the geodesic line. There 
would be no difficulty in performing the construction geometrically; but it would, 
I think, be more convenient to calculate (a) numerically for a given value of V and 
for the successive values of X. Observe that for \’=0 we have (as above) 90°—a=l’, and 
then Tak ==> = pA consequently tan (a) = 5 cot V: but we have also cot 1 =“ cot l, 
so that this equation becomes tan (a) = cotl, or we have 90°—(a)=l; viz. in the 
projection, the geodesic line cuts the equator at an angle =the normal latitude of 
the vertex of the geodesic line. 


The preceding formulæ and results have enabled me to construct a drawing, on 
a large scale, of the stereographic projection of the geodesic lines for the spheroid, 
polar axis =4 equatorial axis. 


Cc. VII. 4 
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